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1 Introduction 

Far-reaching existence theorems, exemplihed by Yau’s solution of the Calabi 
conjecture for Kahler-Einstein metrics, are a distinctive feature of Kahler 
geometry. There are many open problems—the extension of Yau’s results to 
the “positive” case and to Kahler Ricci solitons; existence theory for constant 
scalar curvature and extremal metrics; the study of the associated parabolic 
Ricci flow and Calabi flow equations—making up a very active and challeng¬ 
ing research area today. Alongside this, little attention has been paid to the 
search for explicit numerical solutions: that is, good approximations to the 
metrics treated by the theory. The only steps in this direction seems to be 
the pioneering and recent work of Headrick and Wiseman jH] . In this article 
we discuss another approach to this question, illustrated by some numerical 
results for a particular K3 surface S; the double cover of the projective plane 
branched over the sextic curve x^+y^ + z^ = 0. The Kahler-Einstein metric ut 
on S is characterised by the fact that, after suitable normalisation, the norm 
of the standard holomorphic 2-form 6 is equal to 1 at each point. We 
hnd a metric utg on S, given explicitly in terms of 26 real parameters, such 
that the norm differs from 1 by at most about 1.5% and on average, 
over S, by about .11%. Thus there are good grounds for thinking that cug is 
a very reasonable approximation to the Kahler-Einstein metric, which could 
be applied to investigate any specihc geometric question. The author has not 
yet attempted to investigate such applications in this case, except for some 
discussion of the spectrum of the Laplacian which we which give in Section 
4. We will focus on this one example, although it should become clear that 
the methods could be applied, in a practical fashion, to other cases. The 
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author has some other results for toric surfaces which he hopes to describe 
elsewhere. 

In Section 1 below we describe a general approach to these approximation 
questions for projective algebraic manifolds. This approach is quite different 
from that of Headrick and Wiseman. The distinctive features are 

• The use of metrics furnished by projective embeddings of the manifold. 

• The construction of approximations to the differential-geometric so¬ 
lutions as limits of iterates of maps defined by integration over the 
manifold. 

(It seems likely that these maps, which are dynamical systems on hnite- 
dimensional spaces, can be viewed as discrete approximations to the Ricci 
and Calabi flows.) 

To expand on the the first item above, consider an ample line bundle 
L ^ X over a compact complex manifold X, so for large enough integers k 
the sections of give a projective embedding of X. Suppose G is a positive 
definite Hermitian form on the vector space of sections. Then there 

is a metric h on the line bundle characterised by the fact that if (s^) is 
any orthonormal basis of H^[L^) then the function 

a 

is constant on X. The curvature of the unitary connection on associated to 
h has the form —ikojc where ujg is a Kahler metric on X in the class 27rci(L). 
This metric can also be viewed as the restriction of the standard Fubini-Study 
metric on the projective space (defined by the Hermitian form 

G) to the image of X under the projective embedding. Thus we have a way 
to generate Kahler metrics from simple algebraic data—Hermitian forms on 
hnite dimensional vector spaces. 

The potential utility of this point of view has been advocated over many 
years by Yau, and a fundamental result of Tian na shows that this scheme 
does yield a way to approximate any metric. Tian proved that for any metric 
uj in the class 27rci (L) there is a sequence Gk of Hermitian forms such that 

11^ -^gJI = 0{k~^), 

for any norm on X. Now the dimension of grows like where 

n = dimX, so—from a practical point of view—one might think that good 
approximations by these “algebraic” metrics would require the use of un¬ 
duly large vector spaces (with a Hermitian metric depending on 0{k^"‘) real 
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parameters). This point can be addressed by a refinement of Tian’s result. 
There is a sequence Gk such that 


II^-^gJI =0{,k '') 

for any v. This is explained in the Appendix below. In other words, any 
metric can be very rapidly approximated by algebraic ones. (These approx¬ 
imations are in some ways analogous to the approximations of a smooth 
function defined by truncating the Fourier series.) Now suppose we are in a 
situation where we know (or hope) that a special metric uj (Kahler-Einstein, 
constant scalar curvature, extremal, Kahler-Ricci solution) exists. How can 
we generate a sequence of approximations to uj via Hermitian forms Gk, and 
can we get useful approximations with values of k which are small enough to 
be manageable in practice? These are the questions we take up below. 

The focus of this article is on the explicit numerical results. While there 
is already in place a considerable quantity of rigorous theory to back up these 
methods (and the authors interest in these questions grew out of work on 
abstract existence questions), there are also many points where more theory 
needs to be filled in; but we will not dwell on these here. It is also worth 
mentioning that many of these constructions are closely related to ideas in 
Geometric Quantisation Theory, and in particular the asymptotics of the 
classical limit, but we will not say more about this. 

Lacking any background in numerical methods, the author is not really 
qualified to comment on the comparison of the results here with those of 
Headrick and Wiseman. In one direction, the end-products seem comparable: 
approximations to Kahler-Einstein metrics on very special K3 surfaces with 
large symmetry groups using standard PC’s. The author has not been able 
to directly compare the accuracy of the approximations in the two cases, 
but suspects that the best approximation we present here is comparable 
to the middle range of those achieved by Headrick and Wiseman and that 
by increasing k a bit we would achieve something comparable to their best 
approximation. One fundamental disadvantage of our approach is that it 
is limited to algebraic varieties, whereas Headrick and Wiseman are able to 
vary the Kahler class continuously. One clear advantage, at least in the case 
studied here, is that we only need a few real parameters to store the metrics, 
as opposed to thousands or millions of parameters needed to record the values 
for a lattice approximation. 

The author is grateful to Andres Donaldson for instruction in program¬ 
ming. He would also like to mention the crucial part played by Appendix A 
(“Elementary programming in QBasic”) of the text PP in getting this project 
off the ground. 
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2 General theory 

2.1 Constant scalar curvature 


We return to the discussion of an ample line bundle L ^ X, as in the previous 
section. Thus we consider the relations between two different kinds of metric 
data: 


• Hermitian metrics G on the hnite-dimensional complex vector space 
H\X-L^)- 

• Hermitian metrics h on the line bundle such that the compatible 
unitary connection has curvature —2'KikijJh, where Uh is a Kahler form 
on X. 

As a matter of notation, we write G for the metric on and G~^ 

for the induced metric on the dual space. We also use the notation G = 
(Gap), G~^ = G°‘^ when working with a basis (s^) of 

We have two fundamental constructions. Given a metric h on we let 
Hilb(h.) be the Hermitian metric on 

\\^\\'Hiib{h) = 

where d^h is the standard volume form oj^/nl and R is the ratio 

dimil°(L^) 

^ Vol(X, dfih) 

(which does not depend on the choice of h). In the other direction, given a 
Hermitian metric G on H^{L^) there is a metric FS{G) on characterised 
by the fact that, for any orthonormal basis of we have 

a. 

pointwise on X. The Kahler form uJps{G) is the restriction of the standard 
Fubini-Study metric on projective space (dehned by G) to X —regarded as a 
projective variety. Now we say that a pair (G, h) is “balanced” if G = Hilb(h) 
and h = FS{G). (This terminology was introduced in [0]: the notion had 
been considered before by Zhangjini and Luo [H].) Equally, since in this 
situation either of h and H determine the other, we can speak of a metric h 
on or a metric G on being “balanced”. Let M denote the set of 

hermitian metrics on H^{L^) and dehne T : M —M by 

T(G) = Hilb(F5(G)). 
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Thus, by definition, a balanced metric is a fixed point of T. Now we have, 
from [7j, 

Proposition 1 Suppose that the automorphism group of the pair {X, L) is 
discrete. If a balanced metric G in M exists then, for any point Gq G M, the 
sequence T^{Go) converges to G as r tends to infinity. 

(In the original, preprint, version of [7] this convergence was raised as a 
question. In the version submitted for publication the convergence property 
was stated as a fact, but without a detailed proof. Meanwhile, independently, 
Y. Sano ^1] supplied the detailed argument, in reponse to the original version 

of (Zj.) 

We can spell out more explicitly the definition of the map T. Chang¬ 
ing point of view slightly, let Za be standard homogeneous co-ordinates on 
CP^ and X (Z CP be a projective variety. Start with a positive definite 
Hermitian matrix Gag and form the inverse For ^ G set 

D{z)=J2G'^^ZaZg. 


Then the quotients 

r _ 

~ D{zy 

are homogeneous of degree 0 and can be regarded as functions on CP^, and 
so on X. Then the map T is 

(T(G))^5 = R [ f'ysd/j,, 

where dfi is the volume form induced by the Fubini-Study metric of G. Notice 
that the constant R is chosen so that G'^^T{G).yS = X -|- 1. The geometry 
is unaffected by rescaling the metric G so in practise we work with metrics 
normalised up to scale in some convenient but arbitrary way. 

So far the parameter k has been fixed; we now consider the effect of 
increasing k. The main result of jS] relates balanced metrics to constant 
scalar curvature metrics on X. 

Proposition 2 Suppose that the automorphism group of the pair (X, L) is 
discrete. If X has a constant scalar curvature Kahler metric uj in the class 
27rci(L) then for large enough k there is a unique balanced metric on 
inducing a Kahler metric kujk on X, and 0 ;^ —>■ a; as k ^ 00 . Conversely, if 
there are balanced metrics on for all large k and the sequence ujk converges 
then the limit has constant scalar curvature. 
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Taken together, these two results give a procedure for hnding numerical 
approximations to constant scalar curvature metrics. We choose a sufficiently 
large value of k and then compute the iterate T^{Gq) for some convenient 
initial metric Gq on If k is sufficiently large then the limit as r tends 

to inhnity is a good approximation to the differential geometric solution. 

A crucial ingredient in the proof of Proposition 2 is the Tian-Yau-Zelditch- 
Lu expansion for the “density of states” function. If is an orthonormal 
basis for with respect to the standard metric we set pfc = S 

a function on X which does not depend on the choice of basis. Then if we 
take a fixed metric on L and form the sequence of functions pk with the 
induced metrics on we have 

Pfc ~ r + aik^-^ + 02^-2 + ..., 

where the Oj are local invariants and Oi is l/27r times the scalar curvature. 

2.1.1 A toy example 

We take X to be the Riemann sphere CP^, L to be the line bundle 0(1). 
(Strictly this example does not fit into the framework above, since the pair 
has a continuous automorphism group SL{2,C). However, the theory can 
undoubtedly be extended to relax the condition on the automorphisms, in 
the manner of HDl, SO we will ignore this technicality here.) A basis of 
is given by l,a :,.. .x^ where x is a standard co-ordinate on C. 
We restrict attention to S'^-invariant metrics, for the action x i—>■ on 

the sphere. The invariant metrics are represented by diagonal matrices in 
our basis, specified by the k + 1 diagonal entries, thus 

k 

D = ^ap\x\‘^P. 

p=0 

(Notice that the Op are really the entries of the metric G~^ on the dual space: 
in practice it is easier to work with this rather than G.) The round metric 

on the sphere is given by Op = ^ ^ when H = (1 -|- |xp)^. Due to the 

symmetry we know that for any k the balanced metric is a standard round 
metric on the sphere. 

To make things even simpler, we can consider metrics invariant under 
the inversion x i— *• x~^, so ap = ak-p. Since the geometry is unaffected by a 
overall scaling Op h->• Gap we have just [k/2\ + 1 essential real parameters. In 
the first case, when k = 2, one can evaluate the integrals in the definition of 
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T explicitly using elementary calculus. If we normalise so that Oq = 02 = 1/2 
and write Oi = s the map is represented by the function 


^ ^ scosh ^(s) + — l(s^ — 2) 

2s\/s^ — 1 — 2 cosh“^(s) 

The reader who plots this function will immediately see that the iterates 
T^(so) do indeed converge rapidly to the hxed point s = 1. (The formula as 
written is valid in the range s > 1 but has an obvious continuation to s < 1.) 

We now consider the case k = 6, and evaluate the integrals numerically. 
A typical sequence of iterates is indicated in the next table. 


r 

Oq 

di 

0-2 

03 

0 

.018 

.495 

4.5 

54 

1 

.02833 

.8539 

11.04 

40.16 


.03923 

1.268 

13.38 

34.62 


.05331 

1.645 

14.39 

31.81 

4 

.07150 

1.987 

14.90 

30.09 

10 

.2384 

3.493 

15.57 

25.40 

20 

.7365 

5.400 

15.26 

21.20 

30 

.9488 

5.895 

15.05 

20.21 

40 

.991 

5.983 

15.01 

20.03 

CX) 

1 

6 

15 

20 


This conhrms the convergence that the theory predicts, and similar resuts 
are obtained whatever initial values for oq, oi, 02 , 03 are used. The choice of 
the particular initial values here is made because in this case all the metrics in 
the sequence can be represented as surfaces of revolution in R^, so one obtains 
a vivid representation of the evolution of the geometry of the surface through 
the sequence. The initial values give, roughly speaking, a connected sum of 
spheres joined by a small neck. The first application of T stretches the surface 
into a long “sausage”, which becomes convex after four more application of 
T, and thereafter the sausage slowly shrinks in length to approach a round 
sphere. 

In this example we see that the convergence to the limit is, while steady, 
quite slow. Let ej(r) be the difference between aj(r) and the limiting value 
aj(cx)). From standard general theory we know that there is a constant 
a G ( 0 , 1 ) associated to the problem such that for almost all initial conditions 

ei{r) ~ 

as r —cx). This constant a is just the largest eigenvalue of the deriva¬ 
tive of T at the hxed point, and the vector q is an associated eigenvector. 
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Analysing the sequence, one can see that in this case a is about .8. We also 
see numerically that the corresponding eigenvector has entries approximately 


(1,-2,-1,-4). 


The fact that one gets this eigenvector is easily explained by the SO {3) in¬ 
variance of the problem. The eigenvector corresponds to the second spherical 
harmonic on the sphere. This is given by 



Thus 


and 


/2 


(1- 

\x\ 

2)2 

(1 + 

\x\ 

2)4 

(1 + 

\x\ 

2)6 


(1 - + \x\Y = 1 - 2|x|^ - |a;|^ - A\x\^ - |a;|® - 2|a;|^° + 

A general point to make here is that while the convergence of the se¬ 
quence is slow (since .8 is not much less than 1) it is easy, using this 
standard analysis of the linearisation, to dehne much more rapidly conver¬ 
gent sequences, and of course the same remarks hold in more complicated 
examples. 

This toy model is the only example of the T-iteration which we will con¬ 
sider in this paper. The picture extends readily to the case of toric varieties 
in higher dimensions which, as mentioned before, the author hopes to take 
up elsewhere. 


2.1.2 Extremal metrics 

As we stated above, the results of jH] and [7j are limited, strictly, to the 
case where the pair {X, L) does not have continuous automorphisms, but we 
anticipate that the theory can be extended to remove this condition. In this 
regard, we point out here that there is a straightforward modihcation which 
can be expected to produce numerical approximations to extremal metrics 
in the sense of Calabi |1]. Suppose that {X,L) has a conected Lie group of 
automorphisms and £x a maximal compact subgroup K. Then K acts on 
H^{X, L^) and we restrict to iL-invariant metrics. Then it may happen that 
the sequence T^(Go) does not converge but that there is a sequence gr in 
the complexification of K such that gr(T^{Go)) does converge, as r —>• cx). In 
this case, taking the limits, we again get a sequence of preferred metrics ujk 
which we expect to converge to an extremal metric on X as k ^ oo. While 
the theory here needs to be hlled in, the procedure works effectively in the 
examples that the author has studied numerically. 




2.2 Calabi-Yau metrics 


We now turn to the main topic of this paper: the case when the metric 
one wants to approximate is Ricci-flat. While this case could be treated 
within the constant scalar curvature theory outlined above, there is a slightly 
different, and simpler approach. Let X, L be as before and suppose given a 
hxed smooth volume form du on X. A fundamental result of Yau asserts 
that there is a unique Kahler metric a; on X in the class 27rci(L) which 
realises du as its volume form u^/nl. We propose a method for hnding 
numerical approximations to this metric. The most interesting case is when 
X is a Calabi-Yau manifold, so has a nowhere-vanishing holomorphic n-form 
6 . Then the metric with volume form du = {—i)^6 A 6 is Ricci-flat. 

In a nutshell, we modify the dehnitions of the previous section by using the 
hxed volume form du in place of the Fubini-Study volume djji. Thus we say 
that a Hermitian metric G on iL°(X; L^) is v-balanced if G = Hilb,y(FS'(G)) 
where 

l|s|lLb.(h) = RJ 


The existence of balanced metrics in this context is due to Bourguignon, Li 
and Yau |Sj, see also the recent paper [21. It is closely related to well-known 
results relating moment maps to Geometric Invariant theory, and similar 
extensions were considered by Milson and Zombro HH. In fact suppose that 
u is any positive Radon measure on CP”. We say that a metric Ga^ on 
is u--balanced if 


G^s — R 


Z^Zs 


'CP^ D{z] 


dv. 


When V is derived from a smooth volume form supported on a projective 
subvariety X C CP^ this reproduces the previous dehnition. But we can 
also consider other measures, in particular, sums of point masses (or, in other 
words, the case when subvariety has dimension 0). Let us suppose that the 
measure v satishes one of the following two conditions: 


1. For any non-trivial linear function A on the function log(^^j^) 

on CpTv 

is i/-integrable. 


2. 1 / is a sum of point masses supported on a hnite set Z and for any 
projective subspace P C CP'^ we have 


t/(ZnP) z/(Z) 
dim P -I- 1 X -|- 1 


(The dehnition in the hrst item uses a metric on but is clear that 

the condition does not depend on this choice.) It is easy to see that if v is 
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a smooth volume form on a subvariety X which does not lie in any proper 
projective subspace then the hrst hypothesis holds. 

Now we have 

Proposition 3 If u is a positive Radon measure on CP^ which satisfies 
either condition (1) or (2) above then there is a u-balanced metric on 
and this is unique up to scale. 

We recall the proof briefly. Note that the space of metrics M is the sym¬ 
metric space GL{N + 1, C)/U{N) and there is a standard notion of geodesics 
in M; the images of analytic 1-parameter subgroups in GL{N -|- 1, C). For 
any non-zero vector G we let fiz be the function 

ipziG) = logl^lc-i + ^ logdetG 

on M. The key point is that ipz is convex on all geodesics. Changing by 
a scalar multiple only changes fiz by the addition of a constant. Now given 
our Radon measure v we set, with some abuse of notation, 

This is dehned up to the addition of a constant; for example we can dehne 
^jJz for G CP^ by taking the lift to a vector in of length 1 with 

respect to some chosen reference metric. Now Tj, is also convex, being a 
positive linear combination of convex fnnctions. It is easy to check that a 
metric G is zz-balanced if and only if it is a minimnm of Snch a minimnm 
will exist so long as is a proper fnnction on M, and the convexity yields 
nniqneness. In tnrn, is proper on M if and only if it is proper on each 
geodesic. So what we have to verify is that—nnder either of the hypotheses 
(1), (2)—for each geodesic ray Gt we have cx) as t —cx). Using 

the GL{N + 1, C)-invariance of the problem it snffices to consider a geodesic 
of the form 

Gt = diag(e^“*), 

where I] Aq, = 0 and Aq > ... > Aat. Then 

log (H duz,. 

Consider first case (1). The first coefficient Aq must be positive and we have 
log(^ > log(e^°*|a;on = V + log |a;op. 
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By the hypothesis the term on the right hand side is i/-integrable and so 

> Aotz/(CP^) + Const. 

and we see that —>• cx) as reqnired. In case (2) we write v = Y^ 

for points in CP'^. Let 


a{i) = min{a|a;^*^ ^ 0}. 


Then ~ ct as t —> cxo where 



c = 


Elementary argnments, essentially the same as in usnal Geometric Invariant 
Theory discnssion in ra. ra. show that the condition that c > 0, for all 
geodesics, is eqnivalent to the hypothesis (2). 

Thns we know that, nnder the very mild hypotheses (1) or (2), zz-balanced 
metrics exist. To hnd them, we make the obviuons modihcation to the algo¬ 
rithm of the previous section. We dehne Tj, : M —>■ M in just the same way 
as T but using the measure u. That is, starting with a matrix G, we set 



Notice that this is unaffected by rescaling u. 

Proposition 4 Suppose v satisfies either hypothesis (1) or (2). Then for 
any initial metric Gq the sequence Tf{Go) converges to the v-balanced metric 
as r —> oo. 

The proof is similar to that for the map T, but more elementary. We 
show that decreases the function Tj,; then the conclusion follows from the 
properness of To make the notation simpler we will treat the case of point 
masses, the other case being essentially the same. To prove the inequality 
'^^(T^{G)) < we can without loss of generality suppose that G is the 

metric given by the identity matrix. We can also suppose that the total mass 
Jfi'i is 1. We choose representatives in with = 1. Then 


^^AG) = = 0 . 


We treat the two terms in the dehnition of Tj,: 


>I>„(T(G)) = '£’'‘T^\hc)- + j^logdetr(G), 
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separately. By the concavity of the logarithm function 




Now 

^i/i|^|^(G)-i = ^ PiZ^^hf {T{G)Y^ , 
i i,o:,l3 

but this is 

S -^nGunar^ 

which is 1, since T(G)a/ 3 , T(G)“^ are inverse matrices. So the first term in 
the definition of is less than or equal to 0. For the second term, the 
arithmetic-geometric mean inequality for the eigenvalues gives 


^^^^logdetT(G') < log 


/ Tr(T(G)) \ 

V ^+1 ) 


and the term on the right is zero, since 

i 

Putting the two terms together we have ^y{T{G) < 0, as required. 


To sum up, given a volume form v on our algebraic variety X we have, 
for each /c, an algorithm for hnding the i/-balanced metric on H^{X]L^). 
Moreover, this is robust in the sense that if we approximate v by another 
measure v* which is a sum of point masses—as we have to do in numerical 
integration—the numerical algorithm defined by Ty* will converge to a v*- 
balanced metric provided only that v* satisfies the very mild condition (2) 
(which will happen for any reasonable approximation). Taking the restriction 
of the Fubini-Study metric and scaling by k~^ we get a i/-balanced Kahler 
metric ujk^v on X. Now let k tend to infinity. We expect that the result 
corresponding to Proposition 2 is true, so that the kahler metrics k~^ujk,u 
converge to the kahler metric with volume form u. While this is, again, a 
piece of theory that needs to be hlled in, we will assume it is so for the rest 
of this paper. (In one direction, it is not hard to see that if the k~^ujk,v 
converge the limit must have volume form v. We expect that the harder 
converse can be proved by arguments similar to those in [B].) Thus we have 
another procedure for hnding numerical approximations to Calabi-Yau met¬ 
rics. Experimentally at least, this converges more quickly than the more 
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general “T-algorithm”, which is probably related to the fact that the con¬ 
stant scalar curvature equation is of fourth order in the Kahler potential 
while the Calabi-Yau equation is a second order Monge-Ampere equation. 

To give a toy example, consider again the invariant metrics on CP^ 
and the sections of 0{6). Let u be the volume form of the standard round 
metric. Then, with the same starting point as before, we obtain the iterates 
of Tj, shown in the following table. 


r 

Oq 

di 

02 

03 

0 

.018 

.5 

4.5 

54 

1 

.1395 

2.599 

15.20 

28.12 


.4200 

4.420 

15.54 

23.23 


.6920 

5.297 

15.30 

21.41 

4 

.8568 

5.697 

15.14 

20.61 

10 

.9992 

5.998 

15.00 

20.00 

13 

.9999 

6.000 

15.00 

20.00 


The convergence is much faster. The parameter a governing the asymp¬ 
totic behaviour is now about .42 (with the same eigenvector). 

2.2.1 Refined approximations 

In either variant of the theory, the balanced metrics or //-balanced metrics 
cannot usually be expected to give very close approximations to the desired 
differential-geometric solutions for practical values of k. An analysis of the 
convergence, using the Zelditch expansion, would probably show that the 
convergence is only 0{k~^) or 0{k~‘^). We now return to the issue raised in 
the Introduction of hnding rapidly convergent approximations. We restrict 
the discussion here to the Calabi-Yau case. Thus we suppose that we have 
a metric Go on inducing a Fubini-Study metric ugq on X, and 

that the volume form d/ico = reasonably close, but not extremely 

close to the given volume form u (which we assume to be normalised so that 
the total volumes are equal). We set rj = d/xc/z/; a function on X which is 
close to 1. We would like to “rehne” Gq to get a better approximation. This 
is essentially a linear problem. In the standard differential geometric theory 
we would consider a nearby Kahler metric of the form ugq + idd(j). The 
linearisation of the volume form map is one half the Laplacian of {X,ugo) 
we would take 

0 = -2A-\ri - 1 ) 

where is the Green’s operator. Then, provided that ugq is a good enough 
approxiation for the linearisation to be valid, the metric ugq + i(9(90 would 
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be a better approximation to the desired Calabi-Yau metric and, iterating 
the procedure, we would generate a sequence which converged to that limit. 
Thus the question we address here is how to implement a procedure like this 
numerically, staying within the class of the “algebraic” Kahler metrics. 

Let G be any metric on H^{L^) and Sa be a basis of sections. Then we 
have functions fap = {saiSjs) on X, where ( , ) denotes the pointwise inner 
product on induced by G. Dehne “ 77 - coefficients” 

Val3 = Rj^fapiV - 

(Notice here the diagonal terms rjaa have a particularly simple interpretation, 
when G is the balanced metric and the basis is orthonormal. In that case 
the integral of the positive function f^a is 1/R and the numbers rjaa can be 
viewed as a collection of mean values of rj — l, weighted by Rfaa-) A natural 
criterion to dehne an “optimal” algebraic approximation to the differential- 
geometric solution is that all the coefficients rjap vanish. Thus we say that 
a metric G near to the balanced metric Gq is a “rehned approximation” 
if this occurs. One could hope to prove that, for large enough k, these 
rehned approximations exist and (with a suitable interpretation of “close”) 
are unique. Further, it is reasonable to expect that the resulting sequence 
of rehned approximations would be very rapidly convergent, in the manner 
discussed in the Introduction. 

In more invariant terms, we dehne a vector space Ti.{L^) to be the Her- 
mitian forms on Then, given Gq, we have maps 

^G^{X) 71 : G°^{X) 

dehned in terms of a basis by 

^((^a/?) ^ ( (^apfapj 

and 

(7r(F))„^ = r[ 

Then the condition we are considering is that 7r{f]) G 7i(L^) should vanish, 
which is the same as saying that r] is orthogonal in the sense to the hnite 
dimensional subspace Im(7) C G°°{X). 

The practical question we now face is: how can we hnd the rehned ap¬ 
proximations numerically, starting from the u balanced metrics? As we have 
explained above this is basically a linear problem. If we set G = Go + h then 
we have a map 

V : h {riap{Go + h)) 
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from hermitian matrices h to hermitian matrices rjap- We want to find a zero 
of V and the standard procedure would be to invert the linearisation. This 
linearisation is given by a 4-index tensor Sa/sjs with 


Sa/S'fS fafsif'fS T 2^^f■ys')■ 


If we compute S and invert the corresponding matrix we could dehne an 
iterative procedure which ought to converge to the refined approximation. An 
obstacle to carrying this through is that (even when reduced by symmetry) 
the tensor S has very many components so is relatively hard to compute in 
practice. Thus the author has not yet tried to implement this scheme but has 
used the following simpler procedure instead. Starting with an approximation 
Gq we compute the error matrix E = {rjap) and simply set 


Gi = Gf) — kE 


where k is a suitable positive constant. In fact we compute with the inverse 
metrics Gg ^ and, since E will be small, use the approximation 

G-i = Go^ + kGo^EGo\ 

Iterating this procedure yields a sequence which appears to converges rea¬ 
sonably well, although slowly, to the rehned approximation, see the examples 
and discussion in the next Section. In Section 4 we give some further discus¬ 
sion and theoretical justihcation for this procedure. 

Of course the whole theory sketched here needs to be developed properly, 
and better methods found. But we hope it will yield systematic procedures 
for hnding improved numerical aproximations, starting with the balanced 
metrics. 


2.2.2 Non-zero cosmological constant 

We now consider the problem of hnding approximations to Kahler-Einstein 
metrics with non-zero scalar curvature. Of course this can be htted into the 
constant scalar curvature theory described above, but there is also a natural 
variant of the Calabi-Yau construction above which probably yields a simpler 
approach (although the author has neither attempted to develop the theory 
of this nor studied substantial examples numerically). 

We suppose that our positive line bundle L is either the canonical bundle 
A' of X or its dual K~^. Write p = ±k in the two cases, so our space of sec¬ 
tions is Ar°(X; RP). Given a metric on H^{X] K^), let Sa be an orthornormal 
basis and set 
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4> = Sf). 

This is a section of the bundle ® over X which does not depend on 
the choice of orthonormal basis. Moreover if the sections of generate the 
hbers, which we can suppose is the case, 0 does not vanish on X. Then 
is a well-dehned volume form on X. Using this volume form in place of 
the Fubini-Study form we get the notion of a canonically balanced metric on 
H^{X]KP). Likewise, using this volume form we dehne another variant Tk 
of the map T. Of course we hope that if a Kahler-Einstein metric exists then 
it is the limit of canonically balanced metrics, and that the iterates of 
converge to the canonically balanced metrics. Furthermore, we can combine 
this discussion with that in (2.1.2), in the case when X has continuous au¬ 
tomorphisms, and we can hope to hnd Kahler-Ricci solitons in appropriate 
cases. But we leave all of this as a programme for the future, except to give 
here another toy example. If we take invariant metrics on CP^ and sections 
of 0{6) = K~^ we get the sequence of iterates of T^- 


r 

Oq 

di 

02 

03 

El 

.018 

.5 

4.5 


1 

.0681 

1.714 

14.20 

32.04 


.1860 

3.120 

15.50 

26.39 


.3673 

4.178 

15.57 

23.78 


.5598 

4.900 

15.43 

22.23 


.9821 

5.964 

15.02 

20.07 


.9998 

6.000 

15.00 

20.00 


The convergence is intermediate between the two previous cases, with the 
parameter a about .56. 

3 Study of a K3 surface 

3.1 Geometry 

We now reach the heart of this article: the numerical study of a particular 
K3 surface S. This is the double cover of the plane branched over the sextic 
curve + y^ + = 0. Thus S is dehned by the equation +y^ + , 

where w is a point in the total space of the line bundle 0{3) over CP^. Most 
of the time we work in the affine piece of S which is the subset of C^, with 
co-ordinates {x, y, w), cut out by the equation = x^ + y^ + 1. We hx the 
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nowhere-vanishing holomorphic 2-forni 6^ on S' given in this affine piece by 

^ dxdy dxdy 

w ^/x^ + y^ + 1 

This determines a volume form v = 99. 

The surface S has many symmetries. These are generated by 

• The permutations of x, y, z 

• Multiplication of x, y by sixth roots of unity. 

• The covering involution w i— >■ —w. 

• The antiholomorphic involution given by complex conjugation of all 
coordinates. 


Thus we get a symmetry group T of order 6x6^x2x2 = 864, preserving 
the volume form v. 

We now wish to evaluate the total mass 

Volfs', v) = f V. 

Js 

This is not strictly necessary for our main purpose, but gives a valuable check 
on the accuracy of our numerical calculations in the next section. We can 
evaluate the integral by exploiting the fact that S is an an elliptic surface, 
of a very special kind. To see the elliptic hbration we consider the map 
/ : [x,y,x] I—>■ [x^,y^,z^] from to P^. This maps the sextic branch curve 
to the conic + Z'^ = 0. Thus the covering S' — P^ is the pull back 

by / of the familiar covering of the quadric Q over the plane, branched along 
the conic. Now Q = P^ x P^ is fibred by lines (in two different ways). These 
hbres are the preimages under the covering Q —^ P^ of the lines tangent to 
the conic. The preimages of these lines under / are cubic curves in the plane 
tangent to the sextic branch curve at each intersection point, and these are 
the elliptic curves whose lifts to the double cover yield an elliptic hbration 
of 5. 

To see all of this more explicitly we work in affine coordinates. Let C 
be the plane curve with equation + 1 = 0. Recall that the conic 

-|- /i^ -|- 1 = 0 is parametrised by a rational variable r with 

^ = + = ( 1 ) 

Now for a point (p, q) in the curve C and a complex parameter r set 


yi/3 y ^1/3 ’ 


( 2 ) 


17 



where A = A(r), /x = /i(r) are given by Equation[T] Of course we have to deal 
with the cube roots in Equation|21, so initially we just regard r as varying in 
an appropriate open set O in P^. A few lines of algebra show that if we set 


w = 



then X, y, w satisfy the equation + 1. In other words, we have 

defined a holomorphic isomorphism F from O x O C x C* to an open set 
in S. Further straightforward calculation shows that 


2^/3 

F*{e) = ^ 0 ^ 


where 


dp 


(l_p3)2/3 




dr 

(r(r^ — 1)^/3 


The form 0 is just the standard holomorphic 1-form on the elliptic curve C. 
The form 0 is initially defined only an open set hi in P^. We could introduce 
a covering S —P^ such that 0 lifts to a holomorphic form on S, but there 
is no need to do this because the 2-form 0 A 0 is a well defined integrable 
form on P^. Since we can cover a dense open set in S by taking hi to be a 
cut plane, we see that 

24/3 

Vo\{S,iy) = —IJ, 


where 


/ = / 0 A 0 , j = i ypi 0 A 0. 

Now let W we the wedge-shaped region in C defined by 0 < arg(p) < tt/S. 
A dense open set in the curve C is covered by 18 copies of W (6 rotations in 
C each of which has 3 lifts to C). So 


I =18 [ 0 A 0. 

Jw 

Now it is easy to see that the indehnite integral of (1 — p3^-2/3 j^^aps W to 
an equilateral triangle in C with side length Lj where 


It follows that 


^ 0 dp 

' Jo (l-p3)2/3- 



(3) 


18 








Similarly let Q be the domain in C defined by the conditions |r| < 1,0 < 
arg(r) < 7r/2. This is a fnndamental domain for an action of a gronp of order 
8 generated by r i—>■ and r ^ ir which preserves the form 'ip So 


J = 8 Ip f\ Ip. 

Jq 

The indefinite integral of p) maps Q to another eqnilateral triangle with side 
length Lj say, so 

J = 8—L^j, 

2 


Lj = 


dr 


where 

^ Jo (r(l — 
Finally, the elementary snbstitntion. 


P = 


+ 1 


2/3 


(4) 


shows that 


42/3 

Lj = S^Lj. 

(We can also interpret the snbstitntion of Eqnation 0 as defining a covering 
map from S to F.) Pntting everything together we obtain 

Vol{S, u) = 27Lj, 


where Lj is the one-dimensional integral in Eqnation 01 Evalnating this 
nnmerically one gets Lj = 1.76664, which yields 


Vol(F, u) = 263.000, 

accurate to about 5 or 6 significant figures. 


3.1.1 Linear systems on S. 

We will study the general procedures of Section (2.2) for line bundles 0{k) 
over S: powers of the lift of the hyperplane bundle 0(1) on P^. We will 
consider three cases, when k = 3,6,9. The symmetry group of S acts 
on H^{S] 0{k)) and we can restrict attention to invariant metrics. It is a 
straightforward matter to describe these. There is a natural holomorphic 
section w of 0{3) over S and for /c > 3 we have 

H\S-, 0{k)) = i7°(p2, 0{k)) © 0{k - 3)). 
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The two summands are eigenspaces of the action of the covering involution, 
so must be orthogonal for any invariant metric. The sections of 0{j) over 
are represented, in our affine piece, by linear combinations of monomials 
ioT p + q < j. So we have a standard basis of the space labelled by the integer 
points in a triangle. We will draw this as a right angled triangle but it is better 
to think of it as an equilateral triangle (with a hexagonal lattice) which makes 
the action of the group of permutations of the projective coordinates x, y, z 
apparent. Now considering the action of multiplication by sixth roots of unity 
we see that we can only have a nontrivial inner product between monomials 
x'Py^ and x''y^ if p = r and q = s modulo 6. These inner products must be 
real numbers, due to the symmetry under the antiholomorphic involution of 
5 . 

Now for /c = 3 an invariant metric must be diagonal in our standard ba¬ 
sis. The metric is specified by 4 real parameters a/, a//, anii hi as indicated. 



In other words, the function D associated to the metric is given, in our 
affine coordinates, by 

D = ai\xy\^+aii{\x\^+\x\‘^+\y\^+\y\'^+\x^y\^+\y^x\‘^)+ain{\x\^+\yf+l)+bi\w\^ 


For k = 6 we again have diagonal elements hlling up two triangles, 
invariant under the permutations of x, y, z. This gives us 10 parameters 
a/, an,... ayii, hj, bn, bni according to the scheme indicated: 
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We also have three non-trivial off-diagonal terms corresponding to the 
inner products between These must all be equal, by symmetry, so 

we get one further real parameter C. In other words the matrix of our inner 
product contains a 3 x 3 block 

^ ayii C ^ \ 

C Q-vii C 1 . 

\ C C avii / 

(We could diagonalise this by using a different basis, but to fit in with 
the notation below we will not do so.) 

In other words, the function D is given by 

D = ^ ap^g\x\'^^\y\‘^^ + C{x^lf + x^y^ + + x^ + y^ + y^) + {101“^ ^ bp^q\x\‘^^\y\‘^'^, 

where the coefficients ap^g, bp^g are given by the parameters a/,..., ayii, &/,•••, 
in the manner indicated above. 

Finally we consider the case k = 9. The diagonal elements in the big 
triangle are specihed by twelve parameters aj,, axii in the manner indi¬ 
cated (here we only draw a piece of the big triangle, the remainder follows 
by symmetry). 


bill 
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( 2 / 



am 

an 


avn 

avi 

ay 

aiv 

o^xi 

ax 

O/X 

aviii aviii 


There are now more allowable off-diagonal elements, specified by 6 inde¬ 
pendent parameters Ci,..., Cg. These fit into 3x3 blocks in the matrix of 
the metric of the following forms: 


/ 

avi 

C'2 

c, \ 

/ 

dxil 

c. 

C 4 \ 

/ 

0-VII 

Cg 

Cg \ 


C'2 

ax 

Cs 


c. 

<^IX 

C 5 , 


Ce 

avii 

Cg 

V 

C'l 

Cs 

dxi / 

V 

c. 

Cs 

dix j 

V 

Ce 

Cg 

avii J 


The discussion of the small triangle repeats that in the k = 6 case above: we 
have 8 parameters hj,, byii, C. 

To sum up, we have now specified how invariant metrics on 0{k)) 

are determined by 4,11 and 26 real parameters respectively in the three cases 
k = 3,6, 9. Notice that these are vast reductions on the dimensions of the full 
spaces of metrics dim which are 121, 1444,6889 respectively. 

3.2 Numerical results 

3.2.1 Numerical volume 

The numerical implementation of the algorithm of (2.2) for finding a balanced 
metric is completely specified by the choice of an approximation u* = ^i^zi 
for the given volume form v on S. That is, what we compute is the sequence 
T^*{Ho). Our numerical results are obtained using a family of such approx¬ 
imations which depend on four integer parameters na,, n„, n^: the larger 
these parameters are the more points Zi are used, which should give a better 
approximation, at the expense of extra computation. More details are given 
in the next section, suffice it to say here that nx,np determine the contri¬ 
bution from one chart and nu,nw from another. We should also emphasise 
that our approximation exploits the invariance of the functions we need to 
integrate under the group of symmetries of S. 
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As a first test of our approximations we compute the volume Vol(5', z/q) of 
S in the measure u. We obtain then a family of approximations V{rix, Up, Uu, n^) 
with 

{P'Xi ^pi 'Tlui ^w) ^p) A ^ui)- 

We hnd: 


l/i(10,10) = 265.84 

1 / 2 ( 10 , 10 ) = 6.5557 

l/i(20,20) = 256.61 

1 / 2 ( 14 , 10 ) = 6.3227 

W(30,30) = 256.88 

1 / 2 ( 20 , 20 ) = 6.2575 

Hi (40,40) = 256.70 

1 / 2 ( 24 , 20 ) = 6.2536 


So, for example, 1/(40,40,24,20) = 256.70 + 6.25 = 262.95 and 
V (20, 20,14,10) = 262.93. Recall that the theoretical analysis gave Vol(S', z/) = 
263.000. These hgures suggest that, for values of the parameters similar to 
those above, we can integrate reasonably smooth functions on S with an 
accuracy of about 3 or 4 signihcant hgures. 

3.2.2 The case k = 3 

We now implement our “T,y-algorithm” to hnd the z/-balanced metric for the 
line bundle 0(3) on S. We choose, arbitrarily, the initial metric specihed 
by the four parameters (1,1,1,1). Then we obtain, using the approximating 
measure z/*(20, 20,14,10): 


r 

dj 

an 

am 


m 

1 

1 

1 

1 

1 

12.68 

8.758 

5.257 

2.722 

HI 

13.20 

8.807 

4.987 

2.439 

HI 

13.25 

8.811 

4.959 

2.414 

HI 

13.26 

8.812 

4.956 

2.412 

HI 

13.26 

8.812 

4.956 

2.412 


(Recall that everything is preserved by rescaling the metrics. Here, and in 
the similar results to follow, we normalise the metrics for r > 1, up to scale, 
in an arbitrary way.) Rounding oh to 4 SF, the iteration reaches a hxed 
point after four steps. Taking a hner approximation z/*(30, 30, 20,14) makes 
hardly any change: to 4SF, the hxed point is now 13.26,8.816,4.956,2.415. 
So we have some conhdence that these last values of the parameters give the 
balanced metric on i7°(S'; 0(3)) to high accuracy. 

We now examine the Fubini study volume form p on S' determined by 
this balanced metric. (To hx constants, we take the Kahler form in the 
class 27rci(0(3)).) That is, for given parameters Ua,, n^, n„, we compute 
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the ratio nju at each of the points in the support of u*{nx,np,nu,nuj) (see 
Section (3.4.2) for more details of the calculation). Now we have another test 
of our numerical methods by integrating n/u with respect to the measure 
iy*{nx,np,nu,nuj). By Chern-Weil theory 

[ da = 9.2.4.71^ = 710.61. 

Js 

so the ratio of the total volumes with respect to fi and a is 710.61/263 = 
2.7019. Taking parameters [28, 28, 20,14] we get 

/ -du* = 2.7023 / dv* 

Js V Js 

which is again a fair agreement. Set r] = (2.7019)“^/i/z/, so the mean value of 
7] with respect to z/ is 1 and the deviation of r] from the constant function 1 is 
a measure of the difference between the balanced metric and the Calabi-Yau 
metric on S, normalised by the appropriate scale factor. We compute 

• The maximum value of rj, 

• The minimum value of t], 

• The mean value of |?7 — 1| with respect to a, 

• The distribution function of r] with respect to a. 

More precisely, of course, we compute the maximum and minimum values 
over the support of a z/*(nj;, np, n„, n^), for appropriate parameter values, 
and we compute the mean and distribution function with respect to a*. The 
result (with parameters [28,28,20,14]) is 

max ?7 = 1.496, min 77 = .2501, Mean (|?7 — 1|) = .262 

For the (approximate) distribution function we give in the hrst row of 
the table a collection of ranges and in the second row the percentages of the 
//-volume of S where rj takes values in the given range. 


.25-.4 

.4-. 55 

.55-. 7 

.7-.85 

.85-1 

1-1.15 

1.15-1.3 

1.3-1.45 

1.45- 

3.4 

6.9 

9.6 

11.7 

15.3 

17.3 

15.0 

15.8 

5.1 


We see that jj, is not really a good approximation to 2.7019z/, as we would 
expect with this very low value of k. So we do not proceed to try to hud a 
“rehned approximation” in this case. 


24 























3.2.3 The case k = 6 


We carry through the same procedure as in the previous subsection. With 
parameters [20, 20,14,10] the following table shows the convergence to the 
balanced metric. 


r 

dj 

an 

am 

aiv 

ay 

avi 

avii 

bi 

bii 

bin 

C 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

51.7 

43.3 

33.1 

23.9 

21.8 

15.9 

7.55 

21.8 

14.1 

7.37 

.915 

2 

56.4 

45.9 

33.6 

23.2 

20.8 

14.3 

5.85 

20.3 

12.2 

5.66 

.509 

3 

57.6 

46.5 

33.7 

23.0 

20.6 

13.9 

5.47 

19.9 

11.7 

5.28 

.432 

4 

57.6 

46.6 

33.7 

23.0 

20.6 

13.9 

5.46 

19.9 

11.7 

5.27 

.430 

00 

57.6 

46.6 

33.7 

23.0 

20.6 

13.9 

5.45 

19.9 

11.7 

5.26 

.429 


On the author’s PC each step takea about 5 minutes. Analysis of the conver¬ 
gence suggests that the parameter a is about .22. Increasing the numerical 
integration parameters to [30, 30, 20,16] gives a hxed point, to 4 signihcant 
hgures. 


57.70 

46.61 

33.67 

23.01 

20.56 

13.86 

5.453 

19.95 

11.70 

5.267 

.4324 


We take these parameters as our numerical balanced metric. 

The volume form of the balanced metric yields a function rj with; 

max ?7 = 1.065 , min 77 = .679 , Mean (|?7 — 1|) = .058. 


The distribution function of rj is given by the following table, in which again 
the first row gives the range and the second gives the percentage of the total 
volume lying within the range. 


-.7 

.7-.75 

.75-.8 

.8-.85 

.85-.9 

.9-.95 

.95-1 

1-1.05 

1.05- 

.07 

1.2 

2.2 

3.3 

4.9 

7.9 

14.0 

35.0 

31.4 


This is a much better approximation to the Calabi-Yau metric. The 
function rj is close to 1 over most of the manifold; further exploration shows 
the set where it deviates substantially from 1 is a neighbourhood of the 
branch curve of the double cover, on which rj is small. 

We now move on to search for a rehned approximation, using the algo¬ 
rithm described in (2.2.1) above. Thus we compute the tensor Ti{ri) G 7i(L^) 
as a measure of the size of the error term. If we work with an orthonormal 
basis of sections then the individual matrix entries can be interpreted as 
analogues of Fourier coefficients of 77 . For simplicity we work instead with an 
approximately orthonormal basis given by rescaling the standard monomials. 
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which we justify by the fact that the off-diagonal term C is relatively small 
and we are working close to the balanced metric. Thus for any metric we 


define 


7]{pq) = RapqJ^{7] - 1) 


\q,\2p\y\2q 


r]{pq) = Rbp,qJ^{p - 1 ) 


\x\‘^P\y'^P\w\‘^ 

D 


and 



where R = dimH°(L'^)/Vol(S', z/q) = 38/263. The terms p{pq),p{pq) have 

the same symmetries as the coefficients so we can extract 11 different terms 

Pj, ..., pvii) Vi) ■ ■ ■) Viii) Vc corresponding to the coefficients a/ ..., ayii, bj,..., bjjj, C. 

As described in Section 2, the diagonal terms pi, ■ ■ ■ ,Pni can be interpreted 

approximately as weighted averages of the error p — 1. (The approximations 

we are making here only involve the interpretetation of the data, not the 

actual algorithms.) 

The next two tables display the first 5 steps of the refining procedure 
with the parameter k equal to 2.5. The sequence of metrics is: 


CLj 

an 

am 

aiv 

ay 

ayi 

ayii 

bi 

bn 

bin 

C 

57.70 

46.61 

33.67 

23.01 

20.56 

13.86 

5.452 

19.95 

11.70 

5.267 

.4326 

53.63 

44.15 

33.47 

23.92 

21.62 

15.08 

6.380 

21.73 

13.19 

6.124 

.6059 

51.62 

43.72 

33.51 

24.18 

21.88 

15.26 

6.394 

22.52 

13.55 

6.136 

.6100 

51.01 

43.50 

33.59 

24.34 

22.03 

15.34 

6.347 

23.17 

13.74 

6.096 

.5949 

50.57 

43.35 

33.62 

24.44 

22.13 

15.39 

6.310 

23.61 

13.87 

6.066 

.5785 

50.22 

43.24 

33.74 

24.51 

22.19 

15.43 

6.284 

24.03 

13.95 

6.047 

.5621 


The corresponding sequence of ^^-coefficients, multiplied by 10^: 


H 



Vi 

Vii 

Viii 

Viv 

Vv 

Vvi 

Vvii 

Vi 

Vii 

Viii 

Vc 

-56.1 

-42.9 

-25.2 

-8.0 

-3.6 

10.1 

32.0 

10.7 

25.1 

38.3 

3.0 

-10.7 

-6.9 

-2.6 

1.2 

1.6 

1.7 

-4.8 

11.2 

7.6 

-2.2 

.3 

-5.9 

-3.3 

-.4 

1.3 

1.4 

.8 

-5.1 

9.3 

4.4 

-3.9 

-.2 

-4.4 

-2.3 

.03 

.8 

.9 

.5 

-3.3 

7.7 

2.8 

-2.8 

-.4 

-3.4 

-1.6 

.2 

.5 

.6 

.3 

-1.8 

6.4 

1.7 

-1.9 

-.6 

-2.8 

-1.3 

.4 

.3 

.3 

.2 

-.8 

5.4 

.9 

-1.3 

-.7 



The size of the error terms behaves as follows: 
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r 

Max. 

Min. 

Mean Error(%) 

IT 

1.065 

.679 

5.80% 

1 

1.046 

.804 

2.61% 

2 

1.043 

.830 

2 . 21 % 

3 

1.041 

.843 

1.94% 

4 

1.039 

.853 

1.79% 


1.05 

.860 

1.75% 


We see that for the hrst four steps the error term does decrease, according 
to any of the measures above. The initial rate of decrease is large but this 
soon slows down. At the hfth step the maximum value increases, although 
the mean and minimum improve. The most obvious phenomenon is that the 
coefficient a/ is decreased, along with the “nearby” co-efficients in the big 
triangle, and the co-efficient bj is increased, along with the nearby coefficients 
in the small triangle. This has the effect of increasing the volume form /i near 
the branch curve, and so compensating for the deviation in the balanced 
metric. We repeat the process many times, with various values of k. After a 
while the maximum value decreases again, and for example we achieve at an 
intermediate stage the following metric, which we call cug: 


46.61 

42.69 

34.94 

25.04 

22.56 

15.34 

6.301 

28.95 

13.79 

6.130 

.2854 


This has 


Max Tj = 1.031 , Min. 77 = .898 , Mean (|?7 — 1|) = .017, 
and ? 7 -coefficients, multiplied by 10 ^: 


-.57 

-.11 

.1 

.11 

.04 

-.3 

.25 

.9 

-.2 

.1 

-.2 


The distribution function is: 


-.925 

.925-.94 

.94-.955 

.955-.97 

.97-.985 

.985-1 

1-1.015 

1.015-1.03 

1.03- 

.32 

.45 

1.0 

1.4 

7.9 

38.8 

33.4 

16.7 

.004 


We contine the process still further. The rate of change becomes ex¬ 
tremely slow, and we hnally give up at a metric we will call cug dehned by 
metric coefficients 


43.34 

42.77 

35.84 

25.31 

22.86 

14.88 

6.362 

30.19 

13.78 

6.207 

.1341 


Which has 


Max 7] = 1.050 , Min rj = .909 , Mean (|?7 — 1|) = .0111 
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and ? 7 -coefficients, multiplied by 10 ^: 


-.21 

.011 

.063 

-.026 

.020 

-.070 

-.004 

-.063 

.024 

.005 

-.047 


We take this metric cug as our best estimate for a refined approximation. 
Now there are several points to make. First, the 77 -coefhcients for cug are 
much smaller than those for the balanced metric ujq, by a factor of more 
than 100. However the very slow movement under our primitive algorithm 
means that it is not completely clear that there is a genuine, exact, rehned 
approximation close to cug. Second, the metric cug is probably not significantly 
better as an approximation to the Calabi-Yau metric than the intermediate 
metric cug, and could well be worse. So our overall conclusion is that the 
ideas we developed in ( 2 . 2 . 1 ) are effective, in this case, in generating some 
much better approximations than the balanced metric (such as cug or even the 
fourth step in the process), although a more exact and sophisticated analysis 
of the whole situation is required. 

3.3 The case k = 9. 

The numerical results in this case follow very much the same pattern as for 
k = 6 above, so we will be brief. The crucial difference, of course, is that the 
error terms are smaller. Recall that we now have 26 parameters which we 
display in four-row form as 


CLj 

an 

am 

aiv 

ay 

avi 

O-VII 

aviii 

aix 

ax 

axi 

0-XII 



hi 

bii 

bill 

biv 

bv 

bvi 

bvii 

Cl 

C 2 

Cs 

C'4 

c. 

c. 

0 


We hnd the balanced metric is: 


992.0 

858.9 

682.7 

572.8 

525.7 

401.8 

227.6 

200.8 

176.9 

130.3 

69.46 

19.60 



460.1 

346.2 

218.2 

147.8 

123.8 

68.75 

19.57 

18.45 

19.26 

17.82 

7.275 

8.245 

33.21 

.1903 


The convergence parameter a is about .33. The balanced metric has 

max ?7 = 1.024 , min 77 = .833 , Meand?] — 1|) = .025. 

The largest modulus of the ? 7 -coefhcients is about 34 x 10“^. 

The distribution function is 






-.88 

.88-.9 

.9-.92 

.92-.94 

.94-.96 


.98-1 

1 - 1.02 

1 . 02 - 

1.3 

1.7 

2.1 

3.1 

4.2 

6.6 

12.1 

29.9 

39.1 


Thus over about 70% of the volume, rj lies in the range 1 — 1.025. We 
carry out the rehning procedure as before. The process seems to work better 
than for /c = 6 , in that the error decreases steadily and our best estimate for 
the refined approximation yields the best approximation to the Calabi-Yau 
metric. This approximation, cug is 


875.6 

798.5 

665.8 

584.5 

539.3 

419.5 

239.8 

214.3 

190.4 

140.3 

76.43 

21.76 



540.5 

386.8 

233.3 

161.8 

134.9 

75.53 

21.76 

19.89 

25.89 

19.57 

8.406 

14.44 

31.69 

.5101 


The ? 7 -coefficients (times 10^) are 


-.09 

-.04 

-.001 

-.02 

-.02 

-.03 

-.05 

-.03 

-.04 

-.03 

.005 

-.06 



.16 

-.002 

-.08 

-.02 

-.02 

.01 

-.04 

-.01 

.02 

.03 

.03 

.04 

-.02 

.03 


The error has 

max ?7 = 1.009 , min 77 = .972 , Mean (|?7 — 1|) = .0022, 
and distribution function: 


-.975 

.975-.98 

.98-.985 

.985-.99 

.99-.995 

.995-1 

1-1.005 

1.005- 

.01 

.1 

.3 

.9 

1.5 

46.6 

46.1 

4.5 


So over more than 90% of the manifold rj is within .5% of 1. Notice that 
in passing from /c = 6 to /c = 9, we reduce the mean error in the balanced 
metrics by a factor of about 2, which is roughly (9/6)^. On the other hand 
we reduce the mean error in the rehned approximations by a factor of about 
5. This gives some small support for the idea that the rehned approximations 
should yield rapid convergence, in k, to Calabi-Yau metrics. 

This concludes our account of the main numerical results of this article. 
Notice that the function discussed in the Introduction is essentially 77 “^'^^, 
so the deviation from 1 will be about half that of 77 . 
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3.4 Computational details 

3.4.1 Charts 

To integrate numerically over S we need an appropriate system of coordinate 
charts and choice of good charts is an interesting geometrical question. Sup¬ 
pose 7]) are local holomorphic coordinates for some chart in an atlas, i.e 
with 7]) lying in some bounded domain D C and with the holomorphic 
2-form given hj 6 = F{^,7])d^d7]. Then we want to realise three desireable 
properties: 

• The oscillation of the function F is not too large on D. 

• It is easy to recognize if a point {^,7]) lies in D. 

• The overlaps of this chart with the other charts in the atlas are not too 

small. 

The last property is needed because our integration procedure will involve 
a partition of unity subordinate to the cover, so we need a smooth cut-off 
function Tp^ supported in D and equal to 1 outside the region corresponding 
to the other charts. The contribution from D to the integral of a function / 
on S will have the shape 


[ f7/jD\F\‘^d^d^d7]d7]. 

JD 

If the overlaps are small the derivative of will be large and the numerical 
integration will not be accurate. 

As a first step we restrict to the open set in S' defined by the conditions 

|a;| < 1.2|2:| , \y\ < 1.2|2:|, 

in terms of the homogeneous co-ordinates {x,y,z) on CP^. Clearly S is 
covered by U and its two images under the permutation action, so to integrate 
a T-invariant function it suffices to work in U. Then we can pass to affine 
coordinates, where z = 1 and U is represented as the polydisc |a;|, \y\ < 1.2. 
Now let Ui be the complement in U of the product of the disc \y\ < 1.2 with 
suitable small neighbourhoods of the points where = —1, chosen so that 
the roots of (1 -T x^) can be dehned on Ui. Then 

(i + ( (i + x«)i/« )) ■ 
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w 


Set q = 


Vi+xfi 


so that 


q-l = 


y 


1 + X® 


Now let U 2 C Ui be the region dehned by the condition Re(g) > —1/2.Then 
Ui is covered by U 2 and its image under the covering involution, (which takes 
q to —q). On U 2 we write g = 1 + p® so 


y = p{2 + 


( 5 ) 


where the root (2 is well-dehned on U 2 since we have cut out the 

zero where g = —1 and 2 +p® = 0. The upshot is that we have an open set 
U 2 S on which we can take x and p as complex co-odinates, with y given 
by Equation El and 

tc = (1 + + X®. 

This gives one chart, which we call the “big chart” in S. The holomorphic 
form in these coordinates is 

(2+p6)5/6(l + ^6)l/3'^P'^^- 


Interchanging x and y takes Ui to another open set U[ (Z U say. The union 
UiU U[ covers all of U save for small neighbourhoods of the points where 
X® = — 1 and g® = —1. We want to dehne coordinates on an open set V in 
S covering a neighbourhood of the points x = y = Then U will be 

covered by Ui, U[ and the 36 images of V under the action of multiplication 
by sixth roots of unity. To integrate a T-invariant function it will suffice to 
work in U 2 and V, provided we use suitable invariant cut-off functions and 
take due account of multiple counting by the symmetry group. 

On V we set u = — y^ and take u and w as local coordinates so 


X = 


+ u — 1 


1/6 



This gives our other chart, the “small chart”, on S. in which the holomorphic 
form is 

dudw 

36x^1/^’ 


The tension set up by the requirement discussed above for the two charts can 
now be seen as follows. In the small chart we need |x| not to be too small 
and in the large chart we need 1 -|- x® not to be too small, but the two charts 
must have substantial overlap. We balance these requirements by dehning 


31 










the domain of the big chart to be the region where Re(a;®) > —(.9)® = —.531 
and the domain of the small chart to be where Re(a;®) < —(.7)® = —.117. 
Then in the big chart |1 + a;®| > .47, and in the small chart |a;| > .7. The 
cut-off function has derivative (with respect to x) about 2(.9 — .7)“^ = 10. 

To integrate numerically in the big chart we use hexagonal lattices in 
the complex x and p variables. These allow us to take exact account of the 
residual Z/6 x Z/6 symmetry (that is the integrands are functions of x® and 
p®) saving a factor of nearly 36 in the calculation. We nearly double this 
factor using the further symmetry under complex conjugation. The lattice 
spacings are proportional to where nx,np are integer parameters 

mentioned in (3.2.1). To integrate numeriaclly in the small chart we use 
square lattices in the u and w variables, taking advantage of the residual 
symmetry under complex conjugation and w i—> —w. The lattice spacings are 
proportional to In both charts the precise domains of integration 

are moderately complicated, and so the exact number of points where the 
functions are evaluated to approximate the integrals is not given by a simple 
formula. Writing Ni{nx,np) for the number in the big chart and iV 2 (n„,nu,) 
in the small chart, we have for example 

W(10,10) = 746, iVi(20,20) = 9644, iVi(30,30) = 45481 
N2{10,10) = 9737, iV2(20,20) = 128149, iV2(24,20) = 180435. 

3.4.2 The volume form 

To evaluate the volume form of the Fubini-Study metric dehned by a given 
collection of co-efficients we work in the co-ordinates x,y. Let V+ be the 
complex vector space of dimension [k + l){k + 2)/2 with a basis labelled 
by the integer points in the “big triangle” and R_ the space of dimension 
{k — 2){k — l)/2 corresponding to the small triangle. Our data gives Hermitian 
metrics ( , )± on V± and so a Hermitian metric ( ,) on © I/_. Let 
r{x,y) = r^{x,y) © r_{x,y) be the vector valued function with the entries 
of r+{x,y) equal to x^y^^ and the entries of r_{x,y) equal to wx^y'^, where 
w = y/1 + x® + p®. (Of course w is only dehned up to a sign, but we choose 
a branch of the square root locally.) Write rx,ry for the derivatives of r 
with respect to x, y. Then away from the branch curve tc = 0 the ratio 
V = dpFs/duo of the Fubini-Study form and the hxed volume form 6 /\ 6 is 
given by the formula 


V = 1 1211 ngdet 

/ (r, r) 

{r, Tx) 

{r,ry) \ 

{rx,r) 

{Tx, Tx) 

{rx,ry) 

kc r r 

\ {ry,r) 



{ry,rx) 

{ry,ry) J 
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(6) 



Of course this is the same whichever branch of the square root we take and the 
calculation is completely straightforward. The disadvantage is that we cannot 
use this formula near the branch curve because of the small denominator 
(although in practice the formula seems to be accurate close enough to the 
branch curve for most purposes). A better formula is as follows. Let v± be 
the vector -valued functions, taking values in V± with entries x^y'^. Thus 
r_|_ = and r_ = wv_. Put 


6^ = {v+,wv^),Sy = {v+wv) 


where subscripts denote differentiation. Set 

{Sx,r) 


Scr. = - 


, Sy 6y 




Now write 

/„ = Sx^Jy = 3y^. 


Dehne 



Qx 

||^||2 ((ll'^+ll + (w W )- 

and Qy 

symmetrically. Then 



^ = ^(^1 + C 2 - C 3 + 2 E 4 - 2 E 5 ), 

where 




Cl 

= l«-P(li'5-li’li4lt-l{4.4)lT 


C2 

1^1 2 \\jx:^y Jy^xW ) 


C 3 

= \w\'^\fxQy - fyQx?\ 


C 4 

= Re + \\Sx\\‘^QyW‘^fy) ; 


C 5 

= Re ((4, Sy){f^QyW^ + fyQxw"^)) ■ 


While it appears more complicated, this expression has advantages over 
Equation ini even away from the branch curve. On the author’s PC it takes 
about 1/5 of a second to evaluate this at a given point (when k=9). So to 
evaluate the volume form at 30,000 points—which is the kind of number we 
use for our numerical integration—takes an hour or two. 


33 







4 The Bergman kernel 

Return to the general picture where L —X is an ample line bundle with a 
Hermitian metric and du is a hxed volnme form on X. Then the space of 
sections H^{L^) has an Hermitian inner prodnct. For a point x E X we 
have an evalnation map 

e, : H\L) ^ L,, 

which is represented by the inner product, so there is an element ax E 
H^{L) ® Lx snch that 

s{x) = {s,ax), 

for any section s G H^{L). For x,y E X we dehne 

K{x,y) = \ax{y)\‘^, 
and the associated integral operator 

QK{f){x) = R [ K{x,y)f{y)duy, 

where R = Vol(X, z/)/dim(i7°(L^). (This factor is included to make Qk 
independent of scalings of z/, and Qk{^) = !•) Let be an ortho normal 
basis of H^{L^) with respect to the inner product. Then 


and if we write 


faf3 ("Sa) '5/?) 


(the pointwise inner product over X), we have 


a,(3 


Thus Qk is a hnite-rank operator whose image lies in the hnite-dimensional 
space V C C'^(X) spanned by the fayd-^- the image iiLL^L^) in the notation 
of Section 2.2.1. The restriction of Qk to V gives an endomorphism of V 
with 

Qk{(^'yS h-yS ) R 'y ( f da) faff 

where (,) denotes the inner prodnct on fnnctions. In other words we can 
dehne a endomorphism Q : 7d(L^) 7d(L^) with matrix 


Qal3,'y5 R^f^sfaQi 
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and i o Q = Qj^ o t. One interpretation of Q is that it compares the two 
natnral inner prodncts on V. If we have a given metric on H^{L^) then we 
can identify Ti^L^) with the self-adjoint endomorphisms of and as 

such we have a standard Hilbert-Schmidt norm on 7i(L^) given by 

ll(®a/3)llfl'S = X! l®o/9p- 

On V we have the restriction of the norm and these are related by 

||i(a)||^2 = R{a,Q{a))Hs- 

Notice that, if we regard 7i(L^) as the self-adjoint endomorphisms of 7i(L^) 
then Q has been normalised so that <5(1) = 1. 

To illustrate these ideas, take L —> X to be the line bundle 0(1) over 
CP^ with the standard metric, and standard area form. As in Section (2.1.1) 
we take the usual action on CP^ and restrict attention to the invariant 
part of V, which corresponds to the diagonal matrices in 7i(L^)). With this 
restriction Q is represented hj a {k + 1) x [k + 1) matrix with entries 


Qij 



2k+ 1 ( ‘2k \ 

[^ + J ) 


0 < i,j < k 


(7) 


(see the discussion in (4.3) below). 


4.1 The linearisation of the algorithm 

One way in which the discussion above enters into our theory is in the analysis 
of the linearisation about a balanced metric. Fix an orthornormal basis Sa of 
for the balanced metric. Suppose Gap = 5ap + €ap is another metric. 
Then, to first order in e. 




_ du, 

X) Sp') 


SO we have 

TyiG) = 5aP + €ap + O(e^), 


where 


€ap 7^ / ^ Sp){Sa^ S^')€a^dv^ 

75 


SO e = (5(e). In other words, the linearisation of the map T at the bal¬ 
anced metric is given hy Q. In particular, the largest eigenvalue of Q on the 
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trace-free matrices is the quantity a which determines the asymptotic rate of 
convergence of a sequence T^{Gq), for almost all initial conditions Gq. Thus 
we can estimate this largest eigenvalue, in the examples discussed above, by 
analysing this convergence. For example, on the K3 surface S we estimate, 
by analysing the sequences, that the eigenvalue is approximately .22 when 
k = 6 and .33 when k = 9. 

4.2 Refined approximations and the heat kernel 

Another way in which the operator Q enters our theory is in the algorithm 
we have used for hnding “refined approximations” as discussed in (2.2.1) 
above. Recall that the linearisation of the map which assigns the volume 
form /i^ = cn^/n! to a Kahler metric u is given by one half the Riemannian 
Laplacian, i.e. 

Suppose given any metric G on dehning a metric u on X, and let 

Sa be an orthonormal basis of sections. Consider a small perturbation of G 
to a metric with matrix Sai 3 + in this basis. Then to first order in e the 
induced Fubini-Study metric changes by idd(j) where 

So, to first order, the change in the volume form is 

*^ 0/9 '®/ 3 ))" 

Now given a hxed volume form u on X, write = r]u, where we suppose rj 
is close to 1. As in (2.2.1) dehne 

Vaf3=R (h - l){s^,S/s)du, 

J X 

and consider the variation 

Then the change in the volume form is, to first order, 

'^AQxiv)- 

Thus the algorithm of (2.2.1) will replace an initial error term 77 by a new 
term which is approximately 

^iv) = v-'^^QKiv), 
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so we would like to argue that, for appropriate values of the parameter k and 
with respect to a suitable norm, the linear map hh is a contraction. 

To give evidence for this, we argue that the operator Qk should be related, 
asymptotically as k ^ oo, to the heat kernel on X. Consider the model 
case of a line bundle over C"' with curvature —iu, where u is the standard 
Kahler form (corresponding to the Euclidean metric). Fix a trivialisation of 
the line bundle in which the connection form is ^{J^i^adya — Vadxa) where 
Za = Xa + iya are standard co-ordinates on C”. Then, in this trivialisation, 
the section ao which represents evaluation at 0 is 


so our kernel is 

The Euclidean heat kernel is 


(27r)2- 


H{0,z,t) 



P|2/4t 


so K{0,z) = (27r)“"'Tr(0, z, 1/2). Thus it is reasonable to expect that, on 
a general manifold X the operator Qk will be approximately when 

k is large; so the manifold has very large volume and the local geometry 
approaches the Euclidean case. (Notice that the factor of 27r"^ is accounted 
for by the scaling built into the dehnition of Qk-, since when k is large, by 
Riemann-Roch dimi7°(L^) is approximately (27r)“” times the volume, in the 
metric dehned by the curvature form of L^.) The author has not yet found 
any statement of exactly this kind in the literature but there are results very 
close to this in [H], for example. In any case our present purpose is to give 
a plausible justihcation for the method rather than a rigorous proof. Of 
course, the Laplace operator considered above is that with respect to the 
“large volume” metric, with volume 0{k^). With this discussion in place, we 
argue that near to the Calabi-Yau metric, the operator W is approximately 

IT = 1 - 

2 


Now W is easy to analyse in terms of the spectrum of the Laplacian. On 
an eigenspace belonging to eigenvalue y W acts as (1 — Since 

the function xe~^ has maximum value for positive a:, the operator W 
is a contraction provided that 0 < k < 2e. This is consistent with the 
values of the parameter k found to be effective empirically. These ideas also 
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explain why the “rehning algorithm” takes a long time to get very close to 
the rehned approximation, since the contraction factor for large eigenvalues 
fi is extremely close to 1. 

As a byproduct of these ideas, we can hope to get information about the 
spectrum of the Laplacian of the Calabi-Yau metric from our theory. Let Aq 
be the Laplacian of the metric scaled to have total volume and write 

k' = (dim 

Then we expect that the spectrum of Q approximates that of Thus 

if A is the hrst eigenvalue of Aq we expect that the convergence parameter 
associated to our algorithm is approximately . If, as in this paper, we 

work with L-invariant metrics then we should take A to be the first eigenvalue 
on L-invariant functions. Our estimates = .22, ag = .33 are reasonably 
consistent with this since 

-21og(.22)(38)^/2 = 18.7 

and 

-21og(.33)(83)^/2 ^ 20.2. 

So we expect that A is about 20. (We can also hope to get approximations 
to the eigenfunctions of the Laplacian from the eigenvectors of Q.) 


4.3 Algebraic approximation to the heat kernel. 

We have now explained the importance of the hnite-dimensional linear oper¬ 
ator Q in our theory, and its (probable) relation to the Laplace operator on 
the manifold. Recall that the matrix entries of Q, in terms of an orthonormal 
basis Sa of are 


Qafi^'yS ^ I ('^7’ 


On the face of it, this requires us to evaluate the large number dimiL°(L^)^ 
of integrals over X to hud the matrix. However we can write 


Qcif3,"f5 I ('Sa'5(5) S 

X 


where the products are sections of and ( , ) denotes the hbre metric on 
Let Ti be a basis of 77°(L^^). If we know the integrals 


hj 
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then we can compnte the matrix entries in terms of pnrely algebro-geometric 
data: the prodnct map 




( 8 ) 


Explicitly, if 


then 


^ ' PagiTi 


^ag,'y& ^ ^ PaSiPg'yjI-i 


ij- 




This means that we only need to evalnate approximately 2”'(iV +1)^ integrals 
to hnd the matrix. Moreover these integrals are precisely the integrals which 
dehne the map T for the line bnndle In geometric terms, for any vector 
space V we have the Veronese embedding 


P(V") - P(»"(r)). 


A hermitian metric on V defines a standard induced metric on s^V and, 
up to a scale factor, the Veronese embedding is isometric with respect to 
the Fubini-Study metrics. Thus for X in P(V) we get the same induced 
metric by embedding in P(s^(V)). Starting with the canonical embedding 
in V = we get the canonical embedding in P(if°(L^^))*, which is 

contained as a linear subspace in P(s^V). Thus, starting with a metric G 
on we take the standard induced metric G' on regarded 

as a quotient of the symmetric product. Then the calculation of T{G') is 
equivalent to the calculation of the matrix entries. 

Now suppose it happens that G is the balanced metric for and G' 
is also the balanced metric for This will only be the case in rather 

special circumstances, but for example it holds when X = CP^ with u equal 
to the standard area form. In this case we have T{G') = G' and we can 
hnd the matrix entries purely algebraically, in terms of the product map 
Equation |H1 and the original hermitian metric G. For example for the line 
bundle 0{k) over CP^, restricting to the S'^-invariant metrics, we get the 
matrix entries Qij in Equation [3 above. But in any case we can define 
another endomorphism Q of 7^(L^) by this procedure. That is, we take the 
matrix enties 


Qag 'y PaSiPg'jjl't 


IJ, 




where Itj are the inner products in given by the induced hermitian 

metric, regarded as a quotient of the symmetric square. It is reasonable to 
expect that, when k is sufficiently large, the metric G' is close to the balanced 
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metric and hence that Q is a good approximation to Q. To sum up, starting 
with a hermitian metric G on H^{L^) we have a purely algebraic procedure 
for dehning a self-adjoint endomorphism Q on T-C{L^), and when G is the 
balanced metric (or close to the balanced metric) we can expect that Q is an 
approximation to the heat kernel . 

To illustrate these ideas consider first the case of S'^-invariant metrics on 
3“^. The SU{2) invariance of the problem implies that the eigenspaces of Qk 
correspond to spherical harmonics. Let G [—1,1] be the standard height 
co-ordinate on the sphere and p E 3“^ he the pole where z = 1. The kernel 
function K{p, ) associated to 0{k) is proportional to (1 - 1 - zY and so the 
eigenvalue Xm,k of Q associated to the spherical harmonics of degreee m is 

k 1 

Xm,k 2^+1 J 

where Pm is the Legendre polynomial. It is an exercise in Legendre polyno¬ 
mials to show that 

{k' — 1)... {k' — m) 

^ {k' + 1)... {k' + m) ’ 

where we write k' = k + 1 = dimH^{0{k)), as above. Now set \m,k = 

-2k'\og{xm,k)^ so 

f T \ 

\m,k = -2k' + Jj)} ■ 

From the Taylor expansion of the logarithm we see that 

\m,k = 2m{m -h 1) -h 0{k~Y, 

and the limits 2m{m -|- 1) are the eigenvalues of the Laplacian on the sphere 
of area 27r. For example, we have 

Ai 4 = 4-1- .055, A 2 ,io = 12 -|- .15, A 4 30 = 40 -|- .14. 

This also ties in with observed value of a in Section (2.2), since X 2,6 = 
5/12 = .41666.... 

Finally we consider the balanced metric on the K3 surface 3, with k = 6. 
We restrict to the F-invariant part of We work using a orthonormal 

basis of H^{L^) given by rescaling the standard monomials, apart from the 
triple Here we choose scalars A, B such that A -|- Bx^ + By^ and 

the two similar terms given by permutations are orthonormal. Then we hnd 
that the endomorphism Q on the 11-dimensional F-invariant part of 7i(L^) 
has matrix 
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-.91 -.10 .96 .36 1.19 1.45 \ 

5.71 5.15 2.30 3.13 5.44 2.39 

14.8 13.5 6.46 7.50 13.9 6.56 

10.7 10.2 5.63 5.27 10.5 5.72 

12.7 11.5 6.32 4.34 10.8 6.04 

17.7 16.3 9.62 6.23 15.7 9.35 

16.9 12.5 6.38 16.9 12.4 

14.4 4.01 12.9 14.4 

3.47 7.29 4.35 
19.0 13.7 

14.7 

The matrix is symmetric so we omit the entries below the diagonal. Here 
the hrst basis element corresponds to the off-diagonal term, the next seven 
to the entries in the big triangle and the last three to the small triangle. 

We hnd the first six eigenvalues (ordered by absolute value) of this matrix 
numerically. They are 

1.002, .1956, .05857, .02395, -.002669, .002388. 

The first eigenvalue, 1.002, is a substitute for the exact eigenvalue 1 of the 
matrix Q, and the close agreement is encouraging. The fact that the fifth 
eigenvalue is negative, whereas Q is a positive operator, shows that we can¬ 
not take the approximation this far down the spectrum. For each positive 
eigenvalue y we compute A = —2v^logx and we for the second, third and 
fourth eigenvalues we obtain the A-values 

20.12,34.98,46.00 

respectively. The eigenvalue y = .1956 is in reasonable agreement with 
our previous numerical estimate .22 for the hrst eigenvalue of Q and the 
corresponding estimate 20.12 for the lowest eigenvalue of the Laplacian is 
very close to our previous estimate 20.2 from the observed value of a when 
case k = 9. It is perhaps reasonable to predict, based on this discussion, that 
the next eigenvalue of the Laplacian Aq (on invariant functions) is about 35. 

It would be interesting to test this by repeating the work for k = 9, but the 
author has not yet had time to do so. 

5 Appendix 

Here we discuss the fact stated in the intoduction, that given any Kahler 
metric u in the class ci(L) there is a sequence of “algebraic ” metrics Uk 


( .61 -.33 -.81 -.38 -.84 

3.68 8.38 5.43 4.09 

19.6 13.0 10.6 

8.96 7.63 

9.21 

lO-^x 
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arising from Hermitian forms on with oo^ — oo = o{k^) for all v. 

The proof uses the Tian-Zelditch-Lu expansion and is similar to the ar¬ 
gument in [Hj. We start with Tian’s approximation which, in the notation of 
Section 2, is to take FS o Hilb(a;). (Here we are regarding fc as a parameter 
which is supressed in the notation.) Then 

FS o Hilb(a;) = uj + k~^idd\og{p^) 

where is the density of states function X) kap for an orthonormal basis 
Sq. We know that p^ has an asymptotic expansion 

Puj = 1 -|- fc ^01(0;) -|- k ^Cl 2 (<^) + • • •, 

for certain local invariants a* of the Kahler metric uj. Thus uj — uj = 0{k~‘^) 
and the order k~‘^ term is iddai{uj). Now let a;* = a; — k~‘^iddai{uj) and 
consider the metric FS o Hilb{uj*). Applying the expansion, with smooth 
dependence on parameters, we see that 

F5oHilb(o;*) =uj + 0{k-^). 

We can repeat this process to kill of successively as many terms as we please 
in the asymptotic expansion. The correction terms will become progressively 
more complicated, involving contributions from the derivatives of ai{uj) with 
respect to uj, just as in [0]. In this way we obtain, for any given u, a seqence 
of approximations ujk with 0 ;^ = a; -|- o{k^). A standard “diagonal’ argument 
gives a single sequence with difference o{k'^) for any u. 

There is an elementary argument to prove a somewhat weaker result. 
The construction of a Fubini-Study metric ujh from a Hermitian form H can 
be extended to allow indefinite forms H, so long as H is positive on the 
vectors in corresponding to points of X. It is easy to prove that 

any Kahler metric can be rapidly approximated by algebraic metrics in this 
larger class. It is convenient to assume that L is a very ample line bundle 
over X, so the sections of L give an embedding of X in CP'^. (The argument 
can be extended to avoid this assumption.) For k > 1 consider the standard 
Veronese embedding 

CP^ CP^F 

Let Za denote standard homogeneous co-ordinates on CP'^'' and let I 4 C 
C'°°(CP^) be the vector space of complex-valued functions on CP^ given 
by linear combinations of 

'Z‘a'Z‘b 

\^\2k 

Of course the Za are just the monomials of degree k in the homogeneous 
co-ordinates Za on CP^. 
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Lemma 1 The space Vk is the direct sum of the first k eigenspaces of the 
Laplace operator A^pjv for the standard Fubini-Study metric on CP'^. 


To prove this we take as known the analogous and well-known fact for the 
Laplacians on spheres. The sum of the hrst k eigenspaces for the Laplacian 
on 5'™'“^ C R™ is exactly the space of functions given by restrictions of 
polynomials of degree k on R"^. Now consider the Hopf hbration 
CpA^ This is a Riemannian submersion so eigenfunctions of the Laplacian 
on CP^ lift to S'^-invariant eigenfunctions on the sphere. So the sum of the 
hrst k eigenspaces on CP^ can be identihed with the polyomials in Za-,zF 
which are S'^-invariant. But it is clear that these are just polynomials in the 
products ZaZfi. Separating out the holomorphic and antiholomorphic terms, 
we see that these are exactly the linear combinations of the products ZaZb-, 
as required. 

Now it is a standard fact that if / is a smooth function on a compact 
Riemannian manifold and fk is the projection of / to the sum of the hrst 
k eigenspaces of the Laplacian then fk — f = oiJF) for any z/. Let uo be the 
metric on X given by the restriction of the standard Fubini-Study metric on 
CP^ so 

u = uq + iddf), 


for some smooth function on X. Extend arbitrarily to a smooth function 
on CP^ and take f = e^. Then / is a positive real valued function on CP'^ 
and there is no loss in generality in supposing that the projections fk are also 
positive on CP^. By the Lemma the function fk is a sum 


hABZAZs 

{zl'ik 

and the Fubini-Study metric Uk associated to the form with matrix SAB + hAB 
is a; -f i(9(9(log fk — log /). So Uk — uj is o{k’^) for all z/. 
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